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, The coupling of space-time torsion to the Dirac equation leads to effects on the energy 

levels of atoms which can be tested by Hughes-Drever type experiments. Reanalysis of these 
experiments carried out for testing the anisotropy of mass and anomalous spin couplings can 
lead to the till now tightest constraint on the axial torsion by K < 1.5 • 10~ 15 m — 1 

^ : 1 Introduction 

The geometrical frame for General Relativity is a Riemannian space-time. Within this frame 
one can calculate solar system effects and finds within an accuracy of 10~ 4 that all predictions 
of GR are confirmed by experiment. The equivalence principle which is at the basis of the 
q-( geometrisation of gravity is tested even to much better accuracy (for a review see However, 

on theoretical grounds this geometrical frame may be too narrow, and there are indeed many 
reasons to consider a more general geometrical structure as mathematical description of physical 
space-time. One very prominent generalisation is the Riemann-Cartan geometry which (i) is 
the most natural generalisation of a Riemannian geometry by allowing a non-symmetric metric- 



compatible connection, (ii) treats spin on the same level as mass as it is indicated by the group 
theoretical analysis of the Poincare group, and (hi) arises in most gauge theoretical approaches 
to General Relativity, as e.g. in the Poincare-gauge theory pi, ||, ||] or supergravity ]5[. However, 
till now there is no experimental evidence for torsion. On the other hand, from the lack of effects 
which may be due to torsion one can calculate estimates on the maximal strength of the torsion 
fields. This is the purpose of this Letter: We first calculate that torsion in principle influences the 
experimental outcome of Hughes-Drever type experiments (for a review on these experiments 
see [0). Since no effects were observed we get from the accuracy of these experiments upper 
bounds on the torsion strength. Therefore, by means of a reinterpretation of the Hughes-Drever 
type experiments we obtain the till now most stringent upper bounds on the torsion strength. 

While torsion does not influence the behaviour of macroscopic bodies [M ||] it acts on the 
evolution of spin degrees of freedom and can in principle be measured by determining the 
precession of an elementary spin ||. This spin-torsion interaction also modifies in first order 
of h the trajectory of an elementary particle. Turning around the way of reasoning, it is also 
possible to establish torsion by allowing the spin to behave in a way not predicted by General 



Relativity [1C, 11]. The effect of torsion on spin can also influence the outcome of an interference 
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experiment with neutrons 0, 13]. This fact and the lack of any experimental evidence for a 



coupling has been used in Ref. 13] to pose an upper bound of < 10~ 7 m _1 on the strength of 
torsion fields. Another estimate [14] relates the strength of the coupling of spinning matter to 
torsion to the density of polarised particles. We will show that by means of a reinterpretation 
of Hughes-Drever type experiments we can restrict torsion to < 1CP 15 m _1 . This of course does 
not mean that torsion does not exist or does not play an important role in our understanding 
of gravitation. 

In the following we first derive the non-relativistic limit of the Dirac equation in a Riemann- 
Cartan space-time. The resulting Pauli-equation with additional coupling to the axial torsion 
vector is used to derive the Hamiltonian for the energy levels of a bound two-particle system. 
This determines in the non-relativistic regime the energy levels of a nucleus consisting in a core 
and one valence proton which is the physical system usually taken to perform Hughes-Drever 
type experiments. 

A Riemann-Cartan space-time consists in a metric g^ v (n,v = 0, . . . , 3) and a metric- 
compatible connection Y^ v = {° u } — K pu a where {^} is the usual Christoffel connection and 
is the contorsion tensor related to the torsion S^J 7 through S pu u = = -K^. We 
introduce tetrads /i^ (a, b = 0, . . . , 3) through g pu h p h^ = r]^ with r\ = diag( — h ++). 



2 Dirac equation in Riemann-Cartan space— time 

The Dirac matrices 7^ are defined by the Clifford algebra 7(^7^) = gV v and are connected 
with the standard Minkowski Dirac matrices 7 a fulfilling = rf h by 7^ = h£j a . In Dirac 

representation (m = 1,2,3) 

7 (°) = -if3, 7 m = -i(3a m (1) 



o»\ fl 

a m 0/' P V - 1 
a m \ /O I 



« m = ( I n ), P=[ n 1 ) (2) 



with 



o m J ' 75 = [l 0) ■ (3) 

We use the Dirac equation which is derived from a minimally coupled Dirac Lagrangian [15, 0, M 
(for a recent review see also ||) 

= ih^D^ip + % -K pil p ^$ + map 
{} 

= ihj^D^ip - TiK^^ij) + mcif) , (4) 

where 

= 8^ + I> (5) 
with the spinorial representation of the anholonomic connection 

1 

4" 

= ^D,Khl lbl a -- A K^Kh b vlhl a . (6) 

00.. 

and r „ is the Christoffel part of the covariant derivative and connection, respectively. = 
^e pa up K up a is the axial part of the space-time torsion. 



= tD^ a h b v ibi a 
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3 Non— relativist ic limit: Pauli equation in Riemann-Cartan space- 
time 



In order to carry through a non-relativistic approximation of the Dirac equation coupled to 
metric and torsion we first perform a Newtonian approximation of Riemann-Cartan theory. 
Along the lines of Ilq] we expand the metric with respect to the Newtonian potential in a 
quasi-Newtonian coordinate system (dx° = cdt, i,j = 1,2,3) 

500 = -1 + 2^ (7) 
So* = (8) 
9ij = (l + 2^W (9) 



cr 



The corresponding tetrads are 



h° {0) = 1 + | (10) 



h\ 0) = (11) 
h° m = (12) 

Kn = i 1 -^) 6 ™ (13) 

from which one can calculate the matrices 7^. The Riemannian part of the spinorial represen- 
tation of the anholonomic connection is 

fo = ~^diU (14) 

r* = ~^e i3 k TPd k U. (15) 

The additional axial torsion is taken to be approximately constant at the position of the exper- 
iment. 

We insert the metric, tetrads and connection into the Dirac equation and solve it with respect 
to where we neglect squares of the Newtonian potential and products of U and torsion 

falLy = -ih c (\ - tfdiip - —otdiXJi) + hcK (0) j 5 ^ - hcKi&ij) + (l - 3) 0mc?fy (16) 
C/t \ c J ^2iC \ c J 

In order to perform a Foldy-Wouthuysen transformation of the Hamiltonian (compare |T^1 ) we 
split the Hamiltonian into even and odd parts: 



with 



H = + 8 + (3mc (17) 



O = -ihc (l - a% - ^a%U + hcK (0) -f 5 (18) 
£ = (l - ^) (3mc 2 - hcKiE* (19) 



and get 



I j mcV + ( -|~ A ^ - —K^'iliOrr ) I - I 'Sniy - hrK,^' (20) 
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where we neglected all relativistic corrections and squares and derivatives of the torsion. Pro- 
jection onto the large components and elimination of the rest energy gives the Pauli-equation 
in a Riemann-Cartan space-time 

ih-^-ip = — —Aip - —K (0) aHhdiip - mllib - hcK^ib . (21) 
at 2m m y ' 

(The occurence of c in the last term just means that Ki has the dimension 1/length.) The 
coupling to the electromagnetic field is accomplished by 



ih^-if) = -^—V 2 tp - ed>ib - -^-Hi^ip - mUib - —K^aHfiV^ - hcKi<j l i) = Hib (22) 
ot 2m 2mc m w 



with Vj = di — H^4j- The coupling to rotation Qj can be introduced either by starting with 
the appropriately modified metric components (ffl)and (S) (see [p|]) or by performing a unitary 
transformation of the quantum states obeying ( pl|) to a rotating frame |l9[| , see also [^0| and 
references therein. In each case one gets the additional coupling term £li{L l + \Tmj % ) where L l = 
€ l i k x 3 '(— ihVk) is the angular momentum operator. (Since we are treating the non-relativistic 
limit without second quantisation, we encounter no problems which may arise from a non- 
renormalisability of a Fermi type coupling when inserting into the Dirac equation the field 
equations for torsion within the Einstein-Cartan theory |0|.) 



4 The two— particle— Hamiltonian 

The two-particle Hamiltonian is given by the sum of two Hamiltonians of the above form. 
We restrict from the very beginning to the following case given by the experiment: A nucleus 
consisting in a core with vanishing total angular momentum (J = 0) and a valence proton with 
spin S = ^ and some angular momentum L. The wave function is a function of two position 
variables and an angular momentum variable ip = tpj,Mj(xi, X2,t). Then 

H = H l +H 2 + V (23) 

with 

ex4>-miU (24) 

e 2 4> ~ m 2 U - iJL 2 Hi<? 1 —K (0) a l ihV 2i - TicKio 1 (25) 

m 2 y ' 

where V = V(x 2 — x±) is some binding potential, fi 2 = e 2 h/2m 2 c, and Vi m denotes the U(l)- 
covariant derivative with respect to the coordinate x™. The quantisation axis is defined by the 
external magnetic field Hi. 

We introduce the relative coordinate x and a center-of-mass coordinate X 

x := x 2 — x±, X := x\ H x 2 (26) 

777-1 + m 2 m l + m 2 

and insert the corresponding coordinate transformation into the Hamiltonian ( p3|) . We also split 
the electromagnetic potentials Ai and <p into a part due to external sources Af, (j) e and a part 
due to the charge of the other particle A\, </>': Ai = A\ + A\, (ft = (ff + </>'. We absorb the 
electromagnetic field which is created by the particles itself and which is not connected with a 
derivative into a modified potential V'(x) which depends on the relative coordinate only. Then 



#1 
H 2 



Tr 



2m\ 
2m 2 



V? 



4 



the transformed two-particle Pauli equation reads (in the following we omit the index J, Mj 
characterising the wave function) 



2m 2m rcd c \m 2 m\ J ox 1 



^ , ap/ ^ ap, \ \ dip ih . . e dip 
+— ( e i^(^i) + ^A1(x 2 )) + -A (x)- 



mc v ,v ' lV £ " dX i c 3K 'mdXi 

%Tl ■ f 69 61 \ UiD 

+—A){x) hrr " ^HiU l ip + m 2 £/(x 2 )9? + mil7(xi)^ 

c J \m 2 mi J ox3 

t ^ i • ^ ^ i /" m 2 5 5 \ 

m 2 y ' \m al ! ax 1 / 
+y'(x)(p + ei0 c (xi)v9 + e 2 c (x 2 )^ (27) 

where we defined the total charge e := ei + e 2 , the total mass m := mi + m 2 , the reduced mass 
m red = mim 2 /m, used the Coulomb gauge for the external electromagnetic potential, and ne- 
glected squares of the electromagntic potential and products of torsion with the electromagnetic 
potential. We approximate 

—A?(x 2 ) - —A^xi) ~ —A*(X) + ^■ fiied x k V k A c i {X) (28) 
m 2 mi m rcd n> 

eiA e i (xi)+e 2 A e i (x 2 ) ~ eA*{X) + e Tcd x k V ^(X) (29) 

ei0(xi) + e 2 <\>{x 2 ) « e0(X) + e red ^V^(X) (30) 

miC/(xi) +m 2 C/(x 2 ) £rf mU{X) + ^m red a;VV fe V/£/pf) (31) 
where we introduced the 'reduced charge 



e rcd := m rcd ( — - — ) (32) 
m 2 mi 



and the 'reduced Bohr's magneton' 



/"red := 77— (™>i— + m 2 — ) . (33) 



2mc V m 2 mi 

(In the usually considered case mi — > 00 (very heavy nucleus) we get e rcd — > e 2 and // red — > 
e 2 h/2m 2 c = /i 2 .) If we complete the partial derivatives to covariant derivatives, neglect squares 
of the Maxwell potential, and neglect the internal vector potential A l -{x) which is of the order 
c _1 , then we get 

Hip = -^V 2 xi p--^—A x ip + ^x l V l A c AX)ih dip 



2m 2m rcd x * mc jy ' dXi 

+2^ cd x k V k A c j (X)i^- - toHijtp + V\x)<p 



+e</>*(X)<p + e rcA x l V l 4> c {X) V + mU(X)ip + -m xeA x k x l V k V iU (X)y 

with Vxi = — ^A^(X). This is the final form of the Hamilton operator expressed with 
respect to the relative coordinates x and the center-of-mass coordinates X. 
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Next we extract from this total Hamiltonian that Hamiltonian which describes the energy 
levels of this bound system by freezing the center-of-mass motion and keeping the center-of- 
mass coordinate of the atom at X = 0. In addition, we specialise to cj) e = and gauge away 
constant terms so that we get 

Hip = -^—A x ip + 2^ ed x k V k A e i (0)ih^- i - ^H^ip + V'(x)(p 

1 h 2 ■ d 

+-m red x k x l V k ViU (0)<p - i — K( )a l -^-7ip - TicKia % ip . (35) 

2 777,2 ox 1 

Again we neglected squares of the vector potential. We take a constant external magnetic field: 
A$(X) = \e Uk H l X k . Then the Hamiltonian giving the energy levels is 

h 2 d 

H E = A x - fi red H l 5 i:i e i i k x k ih— - n 2 Hicr l 

zm roc j oxJ 

-—K^aH^j + \m TeA x k x l V k ViU(tt) - ficK i a i + V'(x) . (36) 

777,2 OX 1 I 

If we consider rotation we have to add + \ fi<J l ) where l % is the angular momentum with 
respect to the relative coordinates. The spin-rotation term has been dicussed by Mashhoon 



[19, 20 1 and the Q^P term by Silverman [21|. We get various parts for this Hamiltonian describing 



the energy levels of a bound system: 

He = Ho + H em + -f^Newton + -^torsion (37) 

with 

Ho = --^—A + V'(x) (38) 
2m red 

H cm = -Hi(£jj±r + to<A (39) 

#Ncwton = ^m red x k x l V k ViU(0) (40) 
h 2 ■ d 

^torsion = K (0) aH— - hcKiO 1 (41) 

7772 OX 

For the electric proton-nucleus interaction we have V'(x) = —Ze 2 /x. For the nuclear proton- 
nucleus interaction we have to take some appropriate model for the potential of the nucleus, e.g. 
the harmonic oscillator potential or Wood-Saxon potential. Note that there are no Einsteinian 
effects due to the acceleration VC7. This is in agreement with the equivalence principle: The 
effect of gravitational acceleration can be cancelled by a transformation to a suitable accelerated 
frame and therefore does not influence the energy levels. Ho describes the atom without external 
fields, H em the Zeeman effect. The third Hamiltonian is the usual gravitational interaction with 
the Newtonian part of the Riemannian space-time curvature. The last term describes the 
coupling to torsion under consideration. The first term amounts to a spin-momentum coupling 
which will lead to second order effects only. This is the generalised Pauli-equation for the energy 
levels in a Riemann-Cartan space-time. 



5 Comparison with experiment 

We use the above Hamiltonian to calculate the Zeeman-splitting of energy levels in an atom 
in the presence of torsion. We describe the case which is considered in usual Hughes-Drever 
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type experiments (see e.g. Jfl) namely an atomic nucleus which consists in a J = core and a 
valence proton with angular momentum L = 1. Our quantisation axis for the spin is given by 
the external magnetic field Hi (we can now use /i re d ~ ^2 = ^b)- We are going to calculate the 
shifts in the energy levels due to the interaction Hamiltonian ( |4l] ) describing non-Einsteinian 
effects. We use first order perturbation theory. The unperturbed states \J, Mj) are given by 



|3 3 
I 2 ' 2 



, M Mn ,3 iv ( viii.o>\ U\\l,-l)\ ]3 _ 3) _( 

} ~\ J' \*>*'- \ .fiu i\ r l5 ' ./2| 1)0) I' l2 ' 2 ; - v |i,-i) 



(42) 

The interaction Hamiltonian under consideration ([ll|) has the structure (A\.pi + Ak)cr k . We 
neglect effects due to VjVji7 since these effects give energy shifts smaller than 10 -40 eV which 
is too small to be detectable. We get for the corresponding expectation values 

(H\Ai Pl + A k \l^ = (1,1|A|1,1> (43) 
(H\Ai Pi + A k \l±) = |<1,0|^|1,0> - |<1,1|^|1,1> (44) 

(l-l\Ai Pi + A k \l-±) = hl,-l\A z \l,-l) -hhO\A z \l,0) (45) 



(l.-ll^jtPi + ^fcll,-!) = -(f,-f|A|l,-l) (46) 

where we used that the expectation value for expressions linear in the momentum vanishes. The 
transition frequencies turn out to be 

&*;(§- 1) = |<1,1|^|1 } 1) - |<1 5 0|^|1,0) (47) 

Mh4) = |<i,o|A|i,o)-i(i,i|A|i,i>-i(i,-i|A|i,-i> (48) 

fta,(-i--§) = ^(1,-1|A|1,-1)-^(1,0|A|1,0). (49) 

The matrix elements are (1, l|„4fc|l, 1) = (1, 0|^4fc|l, 0} = (f , — l|^4fc|l, —1} = A z so that we get an 
equal shift Tuo = ^hcK z for all three transition frequencies. The search for such a shift during the 
change of the z-axis with respect to the orthogonal nonrotating Newtonian coordiante system 
amounts to a Hughes-Drever experiment. If space-time torsion will be detected it will lead 



to a diurnal shift of the Zeeman singlet line. However, present experiments (sec [22, 23, [2 4}] ) 
didn't detect any effects. Indeed, the experimental setup of Chupp et al. [23] uses two kinds of 
atoms, 21 Ne and 3 He, where the latter serves as magnetometer standard. Both kinds of atoms 
are subject to the same magnetic field which can be controlled to an accuracy of 5B < 10 -10 G. 
Since both atoms possess different 5-factors, #( 21 Ne) = -0.6619^ B and c/( 3 He) = -2.1276^ 
where /Ub is Bohr's magneton of a nucleon, the Zeeman lines are different. During the experiment 
the energy difference of these two Zeeman frequencies 

£( 3 He) - £( 21 Ne) = ( 5 ( 3 He) - 5 ( 21 Ne)) (fi B B - TicK z ) (50) 

can be recorded. Here K z is a function of the orientation of the quantisation axis. If this axis 
is fixed to the surface of the earth it can exhibit a diurnal time dependence. The accuracy 
of the experiments can be described in terms of an effective variation in the magnetic field B, 
5B < 10 -10 G. Consequently, if the influence of torsion hcK z is larger than fi^dB, then an 
observable effect would occur. If one is going to redo this experiment, searches for the above 
described effect and observes a null-result, then we are lead to the following estimate on torsion 

K z < -z-a^SB < 1.5 • 10" 15 irT 1 . (51) 

he 
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From this new version of an already performed Hughes-Drever like experiment we may get the 
up to now best estimate for the space-components of the axial part of a hypothetical space- 
time torsion. It is not possible to test other parts of the torsion tensor since the Dirac equation 
couples to the axial part only. One needs higher spin equations for a coupling to the trace and 
the traceless part of the torsion tensor. 

Originally Hughes-Drever type experiments are designed to search for possible anisotropics 
of space-time, or, equivalently, anisotropies of the mass of quantum systems. These parts will 
lead to a splitting of the singlet line to a triplet line. Therefore, it is possible to distinguish 
between effects due to mass anisotropy and torsion: while the first cause leads to a splitting of 
the singlet line, the torsion shifts the whole line spectrum in the same way. 



Also the experiments designed to search for an anomalous spin-coupling [25] can be used for 



estimating the strength of the torsion coupling. This experiment has been analysed by Mashhoon 



[20] in order to show that it tests indirectly the spin-rotation coupling. Since space-time torsion 
couples to the spin in the same way as rotation ptf), one also can draw the same conclusions 
regarding torsion. One arrives at similar estimates as above. 

Although Krfy leads to second order effects only in the energy shift, it influences the center- 
of-mass motion via a spin-momentum coupling (compare eqn (|34T) ), Such a coupling can be 
tested with atom beam interferometry. Using a spin flip as described in [19] or in |26|| for testing 



other spin-momentum couplings, we get the phase shift 

5<p = K {0) Al (52) 

where Al is the distance between splitting and recombination of the atomic beam. Note that this 
phase shift is nondispersive and does not depend on the interaction time. If we take an absolute 
accuracy 5<fi < 10~ 2 and Al = 1 m and assume a null experiment, then we get the estimate 
K(o) — 10 _2 m _1 for the time component of the axial torsion. The above phase shift which 
comes from a spin-momentum coupling due to the existence of torsion can be distinguished 
through its mass independence from a similar phase shift due to a violation of local Lorentz 
invariance [E6|. 
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